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ABSTRACT 

From the VIMOS VLT DEEP Survey (VVDS) we select a sample of 16 galaxies with 
spectra which identify them as having recently undergone a strong starburst and subse- 
quent fast quenching of star formation. These post-starburst galaxies lie in the redshift range 
0.5 < z < 1.0 with masses > 1O 9 ' 75 M . They have a number density of 1 x 10~ 4 per Mpc 3 , 
almost two orders of magnitude sparser than the full galaxy population with the same mass 
limit. We compare with simulations to show that the galaxies are consistent with being the de- 
scendants of gas rich major mergers. Starburst mass fractions must be larger than ~ 5 — 10% 
and decay times shorter than ~ 10 s years for post-starburst spectral signatures to be observed 
in the simulations. We find that the presence of black hole feedback does not greatly affect the 
evolution of the simulated merger remnants through the post-starburst phase. The multiwave- 
length spectral energy distributions (SEDs) of the post-starburst galaxies show that 5/16 have 
completely ceased the formation of new stars. These 5 galaxies correspond to a mass flux en- 
tering the red-sequence of pa~+q,psb = 0-0038^oqoi M Q /Mpc 3 /yr, assuming the defining 
spectroscopic features are detectable for 0.35 Gyr. If the galaxies subsequently remain on the 
red sequence, this accounts for 38^^% of the growth rate of the red sequence. Finally, we 
compare our high redshift results with a sample of galaxies with 0.05 < z < 0.1 observed in 
the Sloan Digital Sky Survey (SDSS) and UKIRT Infrared Deep Survey (UKIDSS). We find a 
very strong redshift evolution: the mass density of strong post-starburst galaxies is 230 times 
lower at z ~ 0.07 than at z ~ 0.7. 

Key words: galaxies: high redshift, evolution, stellar content, mass function; methods: sta- 
tistical 



1 INTRODUCTION 

Since iHubblel jl92rj) . the bimodality in the distribution of galaxy 
properties has been one of the great curiosities in the field of as- 
tronomy. Bimodality is observed in galaxy colours, morphology, 
star formation rates and galaxy stellar masses. Modern spectro- 
scopic galaxy surveys have allowed us to quantify the bimodal- 
ity precisely, especi ally with the advent of the Sloan Digital 
Sky Survey (SPSS ) dstrateva et af]|200lHKauffmann et alj|2003bl : 
Baldry et al. 2004). Even at high redshift, evidence f or bimodality 
in the galaxy population i s being sought and f ound (Bund v et all 
120051 : IWillmer et aiteood : iFranzetti et al]|2007» . 

Recent observations have revealed that since a redshift of 
around unity the total mass of stars living in red sequence galax- 
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ies has increased by a factor of two ( iBell et alj|2004) . At the same 
time, the stellar mass density of the blue sequence has remained al- 
most constant. The interpretation is that some blue galaxies migrate 
onto the red sequence after the quenching of their sta r formation, 
whils t the remainder conti nue to form new stars (e.g. iFaber et ail 
1 20071 : lArnouts et alj |2007|) . This quenching of star formation is 
apparently occurri ng in galaxies of i ncreasingly lower masses as 
the Universe ages jBundv et al]|200d) . So-called "dry" (gas poor) 
mergers, with little associated star formation, can increase individ- 
ual galaxy masses w ithin the red-sequence, forming the m assive el- 
lipticals seen today ( Ivan Dokkumll2005l : IBell et al.ll200d) . As well 
as allowing the building of very massive ellipticals, such a sce- 
nario matches the observ ed kinematic and photo metric properties 
of present-day ellipticals jNaab et alj200d,l2007h . 

The overall decrease in global SFR density since z ~ 1 ap- 
pears to be caused by a gradual decline in the mean SFR of galaxies, 
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rather than a decrease in t he number of starb ursts dHammer et alj 
1 19971 ; lNoeskeetai]|2007l ; iTresse et alj|2007t) . However, the pro- 
cess, or processes, responsible for this gradual decline remain to 
be determined. What physical mechanisms cause a galaxy to stop 
forming stars, to turn into a spheroid, and thus to enter the red se- 
quence? There are many competing theories which address one 
or both of these problems. Hot gas stripping (sometimes called 
strangulation) of small galaxies as they fall into overdense en- 
vironments leads to a slow quenching of their star formation 
dLarson et all Il980l : iBalogh & Morrisl |2000| : Ivan den Bosch et all 
120081) . In regions of high external pressure, the cold gas may 
also be stripped from the galaxies by ram pressure stripping 
jGunn & Got3ll972l : ICavatte et a"i1ll994f) . Gas rich major mergers 
may perform both tasks, causing a starburst that is strong enough 
to rapidly exhaust t he supply of fuel, and convert the disk galax- 
ies into spheroid s ( Toomre & T oomrd 1 19721 : iBarnes & He rnquist 
1 19921 ; I Naab & Burkertl l2003h . In this latter scenario, subsequent 
cooling from the inter-galactic medium may cause a disk to reform 
and return the galaxy to the blue sequence. Therefore, the existence 
and efficient coupling to the ISM of additional mechanical energy 
from a central active nucleus (often called AGN feedback) has been 
proposed to prevent refueling and allow the red sequence to buil d 
dSpringel et ai]|2005 atlHopkins et al.ll2007l : lKhalatvan et alj|2008f) . 
Observational evidence for the existence, but above all the relative 
importance of each of these scenarios, remains scarce. 

Much progress has been made on understanding galaxy merg- 
ers since the advent of high resolution computer simulations. Real- 
istic observational comparisons are now becoming easier through 
the application of observational filter functions, point-spread func- 
tions and ra diative transfer o f starlight through dust, to simulation 
outputs (e.g. lLotz et al.|[20o3) . However, many areas of parameter 
space remain relatively unexplored with potentially large implica- 
tions for our overall understanding of galaxy formation. For exam- 
ple, the comparison between simulations of galaxy mergers with 
and without mechanical AGN feedback is rarely focused upon, with 
some groups preferring to focus on starburst feedback, and others 
on the possible impac t of an active nucleus (for an exception, see 
iKhalatvan et~ail2008h . 

In this paper we compare observational results directly with 
a suite of sm oothed-particle hy drodynamic simulations of galaxy 
mergers from ljohansson et"al] d20091) . both with and without AGN 
feedback, by converting the star formation histories of the simu- 
lated galaxies into integrated spectra, upon which the same anal- 
ysis can be performed as upon the observations. Rest-frame opti- 
cal spectra contain a wealth of information about a galaxy's past 
and present star for mation rate, often referred to as the galaxy 
"fossil record" (e.g. IPanter e t al. 2003). In particular, the 4000A 
break and US absorption line strengths constrain the specific star 
formation rate and am plitude of recent bursts of star formation 
dKauffmann et alj 12003 ah . At low redshift, the quality of spec- 
tra and models now allow the analysis of large spectral regions 
dCid Fernandes et alj I200I : iToieiro etail 120071 : Iwild et alj 120071 : 
IKoleva et alj 20081) . At high redshift, spectral indices easily sepa- 
rate bl ue from red sequ ence galaxies without the complication of 
dust dVergani et alj 120081) . At these redshifts, where galaxy spec- 
tra are in general noisier, modern statistical techniques present us 
with the opportunity to vastly increase t he quantity and imp rove the 
quality of derived parameters. Following lWild et al1d2007l) . this pa- 
per presents a principal component analysis of the 4000A break re- 
gion of more than one thousand high redshif t galaxy spectra culle d 
from the VIMOS VLT Deep Survey (VVDS. lLe Fevre et aljboolh . 
We use the information encoded in the spectra to recover the recent 
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Figure 1. The "quality sampling rate" between 0.5 < z < 1.0 for a SNR 
limit of 6.5 and in four bins of 1,4b apparent magnitude. 



star formation history of galaxies at an epoch when the Universe 
was very different from today. 

From an observational perspective, so-called "transition" 
galaxies are of significant interest for pinning down which phys- 
ical processes play the largest role in shaping the evolution of 
galaxies. Classes of galaxies which could be transitioning be- 
twee n the blue and red sequences are "post- starburst" galaxies 
(e.g. iTran et alj|2004 lOuintero et alJl2004l : iLe Borgne et alJl2006l: 



Kavirai et alj|2007h~ "and "green-valley" galaxies (e.g. IMartin et alj 
20071). From the bu ild up of the red sequence in the VVDS, 



Arnouts et al. ||2003) measured the net mass flux which has taken 
place from the blue sequence t o the red sequenc e. This amounts to 
9.8 x 10~ 3 M Q /yr/Mpc ;i for a lChabrier! d2003l) initial mass func- 
tion (IMF), or about 1.4 x 10 4 M Q /yr in the entire VVDS survey 
volume. As we shall show in this paper, galaxies which have un- 
dergone sudden quenching of star formation can be identified in 
VVDS spectra for a period of ~ 0.35 — 0.6 Gyr after the quenching 
event. Therefore, if the entire build up of the red sequence were due 
to a physical process associated with fast quenching, such as gas 
rich major mergers, we could expect to identify a net transitioning 
stellar mass of 6 - 10 x 1O 5 M /Mpc 3 , or about 8.5 - 14 x lO n M 
in the VVDS survey, after accounting for targeting rates. This 
could comprise, for example, a few tens of galaxies of stellar mass 
log(M/M )= 10.5. 

The primary questions that we wish to address in this paper 
are, how much mass do we see entering the red sequence after an 
episode of fast quenching? Are gas-rich major mergers an impor- 
tant mechanism for the build up of the red-sequence since z ~ 1? 
Section 2 describes the VVDS spectroscopic survey and Section 3 
presents our method for identifying post-starburst galaxies. In Sec- 
tion 4 we place our empirical results within a theoretical frame- 
work, through comparisons with simple synthesised stellar popula- 
tions and the star formation histories of galaxy merger simulations. 
The results are presented in Section 5 and compared to the low red- 
shift Universe in Section 6. In Section 7 we discuss the global im- 
portance of post-starburst galaxies for developing a complete pic- 
ture of the evolution of the galaxy population. 

Throughout the paper we assume a cosmology with Ho — 
70kms _1 Mpc~ 1 , O m = 0.3 and Ov = 0.7. 



2 THE VVDS DATASET 

The VIMOS VLT Deep survey (VVDS) is a deep spectroscopic 
redshift survey, targeting objects with apparent magnitudes in the 
range of 17.5 ^ Iab ^ 24. The survey is unique for high redshift 
galaxy surveys in having applied no further colour cuts to minimise 
contamination from stars, yielding a particularly simple selection 
function. In this work we make use of ~ 9000 spectra from the 
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publicly available first epoch public data release of the VVDS- 
0226-04 (VVDS-02h) field w hich covers an area of 1750 square 
arcmins ihe Fevre etail2 005 ). The spectra have a resolution (R) of 
227 and a dispersion of 7.14 A/pixel and a useful observed frame 
wavelength range, for our purposes, of 5500-8500 A. 

The first epoch public data release contains 8981 spectroscop- 
ically observed objects in the VVDS-02h field, 8893 of which are 
within the presently defined survey mask. Of these, 7194 have se- 
cure redshifts (flags 2, 3, 4 and 9, see below), are not classified 
as type 1 AGN and are not secondary objects in the slit. To mea- 
sure each galaxy's star formation history we require coverage of 
the rest-frame wavelength range 3750-4140A, which restricts the 
redshift range of our galaxy sample to 0.5 < z < 1.0. As dis- 
cussed in detail in Section l2. 1 . 3 I below. we impose a signal-to-noise 
(SNR) limit on the spectra greater than 6.5 and magnitude limits of 
19 < Iab < 23, resulting in a final sample of 1246 galaxies. 

In addition to the spectra the dataset includes multi- 
wavelength photometry from diverse sources. Data assembled 
within the context of the WDS survey include deep B, V, R, and I 
photometry from the CFHT/CFH12K camei'43 with limiting mag- 
nitud e s of 26.5, 26.2, 25.9, a nd 25.0, respectively dMcCracken et al.l 
Le Fevre et alj[2004h . Further optical photometry from the 
Canada-France-Hawaii Telescope Legacy Survey (CFHTLSB ex- 
tends into the u and z bands. The field has deep NIR im aging down 
to limits of J;=s 21.50 and Kf=s 20.75 jlovino et alj|2005h and NUV 
and F UV coverage as part of the GALEX mission I Martin et al.l 
l2005h . In order to avoid source confusion due to the large point 
spread function (PSF) of GALEX, UV photometry is based upon 
optical priors (Arnouts et al., in prep.). Infrared photometry in the 
IRAC 3.6, 4.5, 5.8 and 8.0 /xm is available f rom the Spitzer Wide - 
area InfraRed Extragalactic survey (SWIRE, Lonsda le et al. 2003). 
We use the source catalogue as released by the SWIRE tearn_J, with 
detection limits at 5<r of 5.0, 9.0, 43.0, 40.0 /ijy, respectively. In 
this paper the photometric data is used primarily for the purp ose 
of deriving stellar masses, as described in Walcher et alj f2008l see 
Section|2.2|l. 



2.1 Accounting for unobserved galaxies 

In order to calculate absolute quantities from galaxy redshift sur- 
veys, such as galaxy numbers and mass densities, it is necessary to 
account for several effects which cause galaxies not to be included. 
For our work there are three effects: 

• The target sampling rate (TSR): only a fraction of all galax- 
ies in the given magnitude range and survey area are targeted for 
spectroscopic follow-up. 

• The spectroscopic success rate (SSR): a fraction of targeted 
and observed galaxies cannot be successfully assigned redshifts, 
due to either sky line contamination or featureless spectra. 

• The spectroscopic quality success rate (QSR): only a fraction 
of the targeted galaxies with successful redshifts have high enough 
quality spectra from which star formation histories can be mea- 
sured. 

The TSR and SSR are common to many results derived from the 
VV DS, such as lumin osity functions, and are explained in detail in 
e.g. lllbert et al. (2005); only the QSR is new to this work. 



2.1.1 The Target Sampling Rate (TSR) 

As described in lllbert et alj J2005r) . in order to maximise the num- 
ber of targeted objects per VIMOS pointing, the TSR is a func- 
ti on of the size of th e galaxy. We obtain the TSR from Figure 1 
of Ilb ert et aT] J2005l) and weight each galaxy according to its size 
^tsr _ i/TSR(ri), where Ti is the x-radius parameter of each 
galaxy. 



2.1.2 The Spectroscopic Success Rate (SSR) 

The second effect, the SSR, is potentially much more complicated, 
because the success of obtaining a redshift for a galaxy depends not 
only on the apparent magnitude of the galaxy, but also on its red- 
shift and spectral type. The rest-frame spectroscopic features which 
allow accurate redshift determinations depend on the spectral type 
of the galaxy and move through the observed-frame wavelength 
range depending on the redshift of the galaxy. In the VVDS each 
object is assigned a redshift quality flag, which approximates the 
security of the given redshift. Galaxies with flags 2, 3, 4 and 9 are 
deemed securQ. In order to account for the changing SSR as a 
functi on of redshift, an d following the method presented in Section 
3.2 of lllbert et alj (2005J), we estimate t he SSR by makin g use of 
the photometric redshifts calculated by lllbert et all |2006), which 
are publicly available through the Terapix ConsortiurrQ. We calcu- 
late the SSR as a function of redshift and magnitude, as the number 
of objects with secure redshifts (flags 2,3,4,9) divided by the total 
number of objects (with flags 0,1,2,3,4 and 9). Photometric red- 
shifts are used to place the galaxies without secure spectroscopic 
redshifts in the correct redshift bin. 

We calculate the SSR in 4 bins of apparent magnitude 19-21, 
21-22, 22-22.5 and 22.5-23, with bins in Az p i lot of width 0.05. 
Galaxies brighter than 19 and fainter than 23 magnitudes are not 
used in this work, the former because of our lower redshift limit of 
2 = 0.5 and the latter because of the signal-to-noise limitation we 
place on the spectra as described in the following subsection. 

We weight each galaxy in our sample by wf SK — 
l/SSR(z spcc , Iab), where z spec is the spectroscopic redshift of the 
galaxy, and Iab is its J-band apparent magnitude. 

It is important that galaxies with different SEDs are treated 
equally during this weighting procedure. Quiescent galaxies are the 
hardest to assign redshifts to, due to their lack of emission lines. To 
check that the completeness for quiescent galaxies is accurately es- 
timated using the SED-independent method described above, we 
split our sample into high and low SFR subsamples, where SFR 
was measured from SED fitting to multiwavelength photometry 
(see Section f2J\ . We repeated the calculation of wf SR for each 
subsample and found them to be consistent with each other over 
the redshift range 0.5 < z < 1 used in this paper. This redshift 
range is expected to be the least problematic, due to the coverage 
of the 4000A break which is strong in quiescent populations. 



2.1.3 The Quality Success Rate ( QSR ) 

The final selection criterium applied to the VVDS galaxies to create 
the sample studied in this paper is that they have spectra of high 
enough continuum SNR to allow us to extract useful information 



http://cencosw.oamp.fr/ 

2 http://www.cfht.hawaii.edu/Science/CFHTLS 

3 http://swire.ipac.caltech.edu/swire/swire.html 



4 Flag 9: redshift determined from single emission line; Flag 4: 100% con- 
fidence; Flag 3: 95% confidence; Flag 2: 75% confidence filbert et al.l2005l) 

5 ftp://ftpix.iap.fr/pub/CFHTLS-zphot-v 1/ 
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on their star formation histories. In practice, this means a median 
per-pixel-SNFQ larger than 6.5. The SNR of the observed spectrum 
depends on both the apparent magnitude of the galaxy, because all 
galaxies are observed with equal exposure times, and the redshift 
of the galaxy, because the flux of a galaxy varies with wavelength. 

Thus, in a similar manner as for the SSR, we define the quality 
success rate (QSR) to be the total number of galaxies with secure 
redshifts and SNR greater than our limit, divided by the total num- 
ber of galaxies with secure redshifts (flags 2,3,4 and 9). We calcu- 
late this in the same four bins of apparent magnitude used for the 
SSR and in spectroscopic redshift bins of width Az spcc = 0.05. We 
weight each galaxy in our sample by w^ SK = l/QSR(z spcc ,I AB ). 
The w^ 11 is shown in Figure[T]for the four magnitude bins. In Sec- 
tion !5.2l we will show that we recover stellar mass densities consis- 
tent with previous work, despite the addition of this new selection 
criterion. 

Imposing a SNR threshold on the data was found to be neces- 
sary for this work. Using spectra with too low SNR causes consid- 
erable scatter in our spectroscopic indices, which makes difficult 
the robust identification of the post-starburst galaxies in which we 
are primarily interested. This problem is easily identifiable: the in- 
creased scatter in the indices causes the number of post-starburst 
galaxies to increase relative to the total sample, and biases the SNR 
distribution of the candidate post-starbursts towards lower values 
than typically seen in the population as a whole. Therefore the SNR 
limit was chosen empirically, by ensuring that the distribution of 
SNR in our post-starburst galaxy sample was similar to that for 
galaxies of other types, and that our total completeness-corrected 
stellar mass densities for each class of galaxy (i.e. quiescent, star 
forming, post-starburst, starburst) remained constant when the SNR 
limit was increased by a small amount. We found that a per-pixel- 
SNR limit of 6.5 was required for the purposes of our study. 

2.2 Measuring Stellar masses 

A full description of the m ethod used to measure the stellar masses 
of our galaxies is given in Walc her et al. I d2008h . to which we re- 
fer the reader for details. Briefly, we derive the stellar mass M* 
for each galaxy by fitting its multi-band UV-IR spectral energy 
distribution (SEP), to mo del stellar populations created using the 
( iBruzual & Charlotll2003h population synthesis models, assuming 
a IChabrien 1 J2003h initial mass function (IMF). We create a li- 
brary of exponentially declining star formation histories with su- 
perposed stochastic top-hat bursts (see Kauffmann et al. 2003 and 
Salim et al. 2005 for descriptions of similar libraries). The method 
uses Bayesian probability distribution functions as described in Ap- 
pendix A of Kauffmann et al. (2003) to derive physical parameter 
estimates for each observed SED, including errors on the parame- 
ters. For a detailed discussion of the accuracy of stellar m ass esti- 
mates using SED fitting methods see lPozzetti et alj d2007l) . 

2.3 Volume correction 

To account for the changing effective survey volume for galaxies 
of different brig htnesses and SED shapes, we use the V max method 
dSchmidtlf 19681) . This weights galaxies according to the fraction of 
the survey volume in which they could have been observed, given 
their brightness and the shape of their SED. The latter point is 

B Per-pixel SNR is calculated from the median flux over error in the wave- 
length range 5500-8500A. 



important: galaxies with different SED shapes are visible in dif- 
ferent survey volumes and this must be correctly accounted for in 
the V m ax calculation dllbert et ai]|2004) . Using the same library of 
models as used in the stellar mass determination (see above), the 
best fit model was scaled to the observed galaxy magnitude and 
shifted in redshift to determine the total redshift path in which the 
galaxy would have been seen, given the survey magnitude limits. 

The SNR limit alters slightly the effective Iab magnitude lim- 
its of the survey, because very few galaxies with apparent mag- 
nitudes above 23 have spectra with a SNR greater than our limit 
of 6.5. We therefore redefine the survey magnitude range to be 
19 < Iab < 23 and remove galaxies from our sample which 
fall outside of these limits. 

2.4 The survey volume 

The survey area is calculated from the photometric survey mask, 
combined with the boundaries of the spectroscopic survey region. 
The total survey area is 0.472 sq. degrees which corresponds to 
a survey volume in the redshift range 0.5 < z < 1.0 of 1.4 x 
10 6 Mpc 3 . 



3 LOCATING RECENTLY DEAD GALAXIES 

The galaxies in which we are interested are those in which a dom- 
inant A/F star population is visible, implying that the formation of 
O and early B type stars has ceased suddenly. Such galaxies are 
identifiable by their strong Balmer absorption lines compared to 
their mean stellar age as measured by their 4000A break strength. 
The traditional method for identifying post-starburst galaxies in- 
volves the selection of objects with both a strong Balmer absorp- 
tion line (usually H/3 or H<5) to identify a recent starburst, and 
no detectable nebular emission (usually [Oil] or Ha, although 
see Balogh & Morris [2000]) to ensure no ongoing star forma- 
tion. Our selection differs primarily by not placing limits on neb- 
ular emission, as we do not want to bias our sample against nar- 
row line AGN which are enhanced in the post-starburst population 
dYan et alfcOO&IWild et alj200llSch~awinski et al.l2007h . The ad- 
ditional use of the 4000A break strength allows us to select galaxies 
with older and weaker post-st arburst features. 

Following the method o f lWildetalJ ( l2007T) we employ a Prin- 
cipal Component Analysis (PCA) of the spectra around the 4000A 
break to measure the stren gth of the Balmer lines with sufficient 
SNR. In lWildetal.1 J2007b the new spectral indices were defined 
based upon stellar population synthesis models to avoid contami- 
nation of the indices by gaseous emission from the galaxies. Be- 
cause of the low spectral resolution of the VVDS spectra, in this 
work we are unable to accurately mask the Balmer emission lines 
which fall in the center of the Balmer absorption lines in which we 
are interested. Therefore, we choose to perform the PCA directly 
on the VVDS spectra themselves. In Section|4]we explain how we 
compare a posteriori to model stellar populations created with a 
population synthesis code. 

3.1 Calculating the spectroscopic indices 

All galaxy spectra in our sample are corrected for Galactic ex- 
tincti on assuming a uniform E(B— V)= 0.027 dMcCracken et al.l 
2003), moved to the galaxy rest-frame and interpolated onto a com- 
mon wavelength grid. Each spectrum is normalised, by dividing 
by the median flux between 3750 and 4140 A. The mean spectrum 
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Figure 2. The eigenspectra of the WDS galaxy sample in the rest-frame 
wavelength range 3750-4140 A. The spectral indices used in this paper refer 
to the amount of each eigenspectrum present in an individual galaxy spec- 
trum. The first index is related to the shape of the spectrum and equivalent 
to the common D n 4000 index. The second index contains the Balmer ab- 
sorption line series from H<5 to H9, together with the distinctive decrease in 
flux bluewards of 3850 A seen in late B to F-star spectra. 



is calculated and subtracted and a PCA is then performed on the 
residuals. Because PCA is a least-squares process, eigenspectra cal- 
culated directly from data can easily be dominated by outliers in 
the input sample of spectra. In order to obtain robust eigenspec- 
tra that are representative of the majority of the population, w e 
use an iterative and robust PCA algorithm teudavari et al.ll2008h . 
The algorithm essentially replaces the least-squares minimisation 
solved by traditional PCA with the minimisation of a new robust 
function of the data, which employs a robust Cauchy-type func- 
tion to limit the impact of outliers. Th e eigenspec t ra ar e presented 
in Figure [2] As described in detail in Iwild et all J2007t) . the first 
eigenspectrum relates to the strength of the 4000A break, and is 
equivalent to the well-known index D n 4000. This has been shown 
to correlate most strongly with the specific star formation rate 
of the galaxy (SSFR, star formation rate divided by stellar mass, 
iBrinchmann et alj|2004h . The second index contains the Balmer 
absorption lines, together with the correlated decrease in flux blue- 
wards of 3850 A found in late B- to F-type stars. 

We calculate the principal component amplitudes (PCs) of all 
galaxies in our sample by projecting the spectra onto the eigenspec- 
tr a. The projection is perform ed using the "gappy-PCA" procedure 
of lConn ollv & Szalav ( 1999), which weights pixels by their errors 
during the projection and gaps in the spectra due to bad pixels are 
given zero weight. We additionally allow the normalisation of the 
spectra to vary as a free p arameter in the pro jection (G. Lemson pri- 
vate communication, see IWild et al.ll2007h . The PCs represent the 
amount of each eigenspectrum present in a galaxy spectrum and 
are our new star formation history indices. Errors are calculated 



during the projection of the spectra onto the eigenspectra, using 
the error arrays of the individual WDS spectra. These errors are 
purely statistical, however, in Wild et al. (2007) we show that for 
SDSS spectra they compare well to the scatter observed between 
duplicate observations. 



3.2 Defining the post-starburst galaxy sample 

In Figure|3]we show the distribution of PC 1 and PC2 for our sample 
of WDS galaxies. To guide the eye of the reader, and to aid in fur- 
ther discussion, we have split the sample into four "classes", with 
boundaries empirically defined from the distribution of points. We 
would like to emphasise that in reality the star formation histories 
of galaxies form a continuous distribution, and any separation into 
discrete classes is an over simplification, greatly reducing the infor- 
mation contained within the full distribution. The primary division 
of our sample is into "quiescent" and "star-forming" galaxies on the 
right and left, based on the value of PCI (i.e. the strength of their 
4000A break). The precise positioning of the boundary is arbitrary, 
and care should be taken when comparing our results for "red" and 
"blue" galaxies to those derived from different observations. To the 
bottom left, systems with very blue continua and very strong emis- 
sion lines are found, which we class as "starburst" galaxies. 

Finally, to the top center we find the "post-starburst" galax- 
ies, with stronger Balmer absorption lines than expected for their 
4000A break strength. We define "post-starburst" galaxies to have 
-0.6 < PCI < 1.0 and PC2> 0.6. While these boundaries were 
defined empirically from the distribution of data (this class are clear 
outliers in PC2), we will later compare with model stellar popula- 
tions to justify their classification as "post-starburst". We find that 
18 galaxies lie in this region at > la confidence. The spectra of 
these galaxies are shown in Figure [4] The selection of these 18 
galaxies is robust to changes in the exact method used to create the 
eigenspectra. 



3.3 Estimating completeness limits 

It is important to estimate the completeness of the survey in terms 
of galaxy stellar mass. This of course depends upon the mass-to- 
light ratio and SED shape of the galaxies in question: flux limited 
surveys are generally complete to lower mass limits for star form- 
ing galaxies than for quiescent galaxies. While the completeness in 
a particular photometric band can be judged empirically, this is not 
possible for the completeness in derived parameters, such as stellar 
mass and star formation rate. We therefore use an alternative ap- 
proach to estimate the mass completeness of our sample, described 
in detail in W08. 

Because we are primarily interested in the post-starburst pop- 
ulation, with mass-to-light ratios between those of star-forming and 
quiescent galaxies, we derive the mass limits relevant for this pop- 
ulation alone. From the model dataset described in Section I2T2I we 
select all model galaxies which have undergone a starburst in the 
last Gyr and lie within our PSB region in PC 1/PC2. For a given stel- 
lar mass, we determine the fraction of model post-starburst galaxies 
that lie within our survey magnitude limits. We find that 50% of the 
model post-starburst galaxies at z = 0.6 with log(M/M Q )> 9.5 
would be detected, and log(M/M Q )> 9.9 atz = 0.9. The stochas- 
tic model library assumed in this analysis has a higher fraction of 
strong post-starburst galaxies than appear in the real data, i.e. while 
the observational PC1/PC2 space is fully covered by the models, 
the relative space density of the different galaxy types is not reli- 
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Figure 3. Left: The first two principal components (PCs) for all galaxies in our VVDS sample. PCI is equivalent to the well known index D n 4000, PC2 
represents the excess (or lack) of Balmer absorption. For analysis purposes, the sample has been split into quiescent (red), starforming (cyan), star-bursting 
(blue) and post-starburst (orange) classes. Post-starburst galaxies are defined to lie within the dotted box indicated by at least la. Those with SSFR< 10 /yr 
are circled. The median errors of the whole sample (black) and the post-starburst galaxies alone (orange) are shown in the top left. Right: The volume density 
of galaxies with log(M/M„)> 9.75 in 0.5 X 0.3 bins inPCl/PC2. 
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Figure 4. The spectra of post-starburst galaxies culled from the VVDS gal 
PCA fits to the spectra using 10 components. 

able. We therefore confirmed that the mass completeness limits did 
not vary greatly with burst mass fraction or burst age. 
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sample in the rest wavelength range 3750^4140A. Overplotted in red are the 
4 COMPARISON WITH MODELS 



In the following sections, all results will be quoted for a mass 
limit of log(M/M Q )> 9.75. We have confirmed that increasing 
this limit to log(M/M Q )> 10.0 does not alter our conclusions. 
Figure [5] (right) shows the number density of galaxies in bins of 
PCI/2 for galaxies with log(M/M Q )> 9.75, corrected for survey 
volume and sampling effects. 



In order to understand the particular properties of a starburst which 
would lead to a post-starburst galaxy in our sample, and to mea- 
sure the timescale during which a post-starburst galaxy is visible, 
it is necessary to compare our data with model stellar populations. 
We begin by creating simple toy model starbursts and subsequently 
turn to the analysis of smoothed particle hydrodynamic (SPH) sim- 
ulations of merging galaxies. 
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Figure 5. By projecting evolving model stellar populations onto the WDS eigenspectra, we can measure quantitative properties relevant to our PSB galaxies. 
After continuously forming stars for 5Gyr, the toy-model galaxies undergo a starburst which declines exponentially with decay time r^cay Left: The 
maximum PC2 attained by the models after the starburst has occurred, i.e. a measure of the strength of the PSB features. Different lines represent models with 
bursts of different T c i ecay . Right: The time spent in the PSB phase as a function of burst mass fraction. 



4.1 Toy model galaxy tracks 

For a galaxy to reveal strong B aimer absorption lines, star forma- 
tion must switch off rapidly resulting in an excess of longer lived A- 
and F-stars over the hotter O- and early B-stars. This does not nec- 
essarily have to be preceded by a short-lived starburst, the sudden 
truncation of ordinary star formation will suffice, although the sig- 
nature will of course be weaker. The two parameters of importance 
for the strength and longevity of the post-starburst phase are the 
mass fraction of stars formed in a starburst and the decay timescale 
of the starburst. In this terminology, a truncation model then simply 
has a mass fraction of zero per cent. 

With this in mind we have created a library of model galax- 
ies which undergo exponentially decaying bursts, occurring after 
5 Gyr of continuous star formation. These bursts are parameterised 
by their decay time and burst mass fraction, i.e. the mass of stars 
formed while the star formation rate is above the continuum level. 

Converting the star formation histories into spectra suitable 
for comparison with the WDS dataset requires several steps. 
We input the star format ion histories into the GALAXEV code 
jBruzual & Charlotl l2003). using Chariot & Bruzual (in pr epara- 
tion) simple-stellar-populations (SSPs), a IChabrieil J2003h IMF 
and assuming solar metallicity. The evolution of the model galax- 
ies in PC 1/2 is not gr eatly affected by the change between 
iBruzual & Chariot! £2003) and Chariot & Bruzual (in preparation) 
SSPs, which involves both a change in stellar libraries and in the 
stellar evolution tracks. Solar metallicity was chosen as being suit- 
able for our post-starburst galaxy sample, given their stellar masses 
(Section l5.1l >. but the evolutionary tracks are insensitive to moder- 
at e variations in m etallicity. The two component dust prescription 
of ICharlot & Fall d2000h is applied to the continuum light as im- 
plemented in the GALAXEV code. Again, the inclusion or not of 
continuum dust has little effect on the time spent on the PSB se- 
quence, or the maximum PC2 observed. This is because the domi- 
nant A/F star population in the post-starburst spectra is assumed to 
already have emerged from the dense stellar birth clouds. 

Next, because the resolution of the WDS spectra is too low to 
mask emission lines in the data, we must add emission lines to the 
model spectra. We convert the rate of ionising photons predicted 
by the stellar population synthesis model {Qi 0n ) into Ha line flux 
using: 

F Ha [erg/s] = 0.45 * E Ha [erg] * Q io „[s _1 ] (1) 



where Eho is the energy of a Ha photon and 0.45 is the fraction 
of ionisations which lead to the emission of an Ha recombina- 
tion photon (case B, Te = 1 x 10 4 K and Ne = 1 x 10 4 cm~ 3 , 
IStorev & Hummed 19951) . Th e remaining B aime r line fluxes are de- 
rived from the ratios given in lOsterbrockl Jl989). The stellar popu- 
lation synthesis model includes the post-AGB phase of stellar evo- 
lution which is believed to be imp ortant for producing ionising pho- 
tons in older stellar populations (Bine tte et al.ll 19941) . We note that 
no contribution to the emission lines from obscured, narrow-line 
AGN is included in the models. Any substantial infilling of the 
Balmer absorption lines in this way will prevent us from identify- 
ing post-starburst galaxies. AGN with strong enough lines to have 
a significant effect are extremely rare however, and their absence 
from our sample will not affect our final results. 

Just as dust attenuation is applied to the model stellar con- 
tinua, the li nes are similarly at tenuated using the attenuation law as 
described in lWild et al] J2007t) : 

This law is i nspire d by the two-component dust model of 
ICharlot & F all ( 2000). For the birth clouds, from which nebular 
emission is believed to originate and which contain a fraction 1 — fi 
of the dust, these authors adopted for simplicity an attenuation law 
to match that of the diffuse interstellar medium (i.e. of the form 
A -0 ' 7 ). Our slightly steeper curve, resulting in slightly increased at- 
tenuation, is motivated by the fact that dense clouds around young 
stars have a more shell-like geom etry than the diffuse interstellar 
medium (see Ida Cunha et al.ll2008l for more details). The attenua- 
tion of the model emission lines by a moderate amount ( Av = 1.0) 
alters slightly the evolutionary tracks in PC1/PC2 of the model dur- 
ing the starburst phase. 

Our ability to detect a post-starburst galaxy depends on the 
maximum PC2 reached during the evolution of the stellar popula- 
tion, and the time spent in the post-starburst phase (/psb)- These 
factors depend in turn on the strength of the proceeding starburst, 
together with the timescale of the subsequent decay of the star 
formation rate (rdccay)- Figure [3] shows that the maximum PC2 
reached by the WDS galaxies is around 1.0, with 5/18 having 
PC2 >0.9. Using our toy models, the left panel of Figure [5] shows 
the maximum value of PC2 attained by a model galaxy after the 
burst has occurred, as a function of the burst mass and burst decay 
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time. The models show that the descendants of starbursts with de- 
cay times < 0.1 Gyr and burst fractions above a few percent will 
appear clearly in the post-starburst phase. The right panel of Figure 
[5]shows the time spent in the post-starburst phase (tpsp,)- This time 
is short, <0.15 Gyr, for burst mass fractions < 5%, and reaches a 
maximum independent of decay time of ~0.6 Gyr for mass frac- 
tions of >30%. We will use this value as an upper limit on tpsB in 
later sections. 

We conclude from this section that the majority of the pro- 
genitors of our PSB galaxies underwent a recent strong starburst, 
followed by a rapid truncation of the star formation. A small frac- 
tion may however be simply the result of rapid truncation of their 
ongoing star formation. 



4.2 Simulation galaxy tracks 

The next step in our comparison wit h models involves a s ample 
of 79 merger simulations presented in Johansson et al. ( 2009j, here- 
after JNB09). 



4.2.1 Simulation details 

The simulations were perf ormed using t he entropy conserving 
TreeSPH-code GADGET-2 dSpringelll2005t) . which includes radia- 
tive cooling for a primordial mixture of hydrogen and helium to- 
gether with a spatially uniform time-independent local UV back- 
ground. Star formation and the associated supernova feedback is 
impleme nted using the sub-re s olution multiphase model devel- 
oped by Springel & Hernquist (2003). We model the feedback 
from black holes foll owing the effective subresolution model of 
ISpringel et al. (2005b), in which the unresolved accretion onto the 
black hole (BH) is related to the resolved gas distribution usin g 
a Bondi-Hoyle-Lyttleton parameterization (Bondi & Hovle 19441) . 
Further details concerning the feedback implementations and pa- 
rameter choices can be found in JNB09. In addition, for the pur- 
poses of this paper a new suite of simulations was performed in 
which the BH feedback was switched off. Such a comparison is 
important because the BH feedback can cause star formation to 
completely shut down after the merger, which may plausibly af- 
fect both the time the merger remnant spends in the PSB phase, 
and the strength of the post-starburst features. 

The simulations include mergers of gas-rich disks (Sp-Sp), of 
early-type galaxies and disks (E-Sp, mixed mergers), and mergers 
of early- type galaxies (E-E, dry mergers). The progenitor galax- 
ies have a range of virial velocities, « v ir of 80, 160, 320 and 
500km/s, which determines their masses and sizes (see Eqs. 1 & 
2 in JNB09). All model galaxies contain an exponential disk (d) 
component, together with a s tellar bulge (6) e mbedded in a dark 
matter halo modelled with the Hernquistl dl990h profile. The stellar 
bulge has a fixed mass ratio of M£ — 0.01367M v i r = l/3Md,tot, 
where the total disk mass is a sum of the disk stellar and gas mass 
M d ,tot = + M d , g as = 0.041M v ir, the disk stellar mass 
is = 0.041(1 — /gas)M v ir, and M v j r is the virial mass of 
the galaxy. The initial gas fraction / gas takes values of 0.2, 0.4 
and 0.8. The galaxies approach each other on parabolic orbits, 
which are generally motivated b y statistics from N-body simula- 
tions ( Khochfar & Burked 120061) . The resulting merger remnants 
have total stellar masses between log(M/M )= 9.5 — 12.5. Both 
major (1:1) and minor (3:1) merg ers were simulated w ith varying 
orbital and initial disk geometries iNaab & Burkertl2003l) . All sim- 
ulations were evolved for a total of t = 3 Gyr, with the merger 



taking place after approximately 1.5 Gyr, using the local Altix 3700 
Bx2 machine hosted at the University Observatory in Munich. In 
this paper we will primarily concentrate on mergers with prograde 
in-plane orbits ("GO" , i.e. the galaxies collide directly edge on 
iNaab & Bur kert 2003). We will additionally show an example in 
which galaxies rotate in a retrograde sense, with the disks tilted by 
30 degrees with respect to one another ("G7"). We will focus ex- 
clusively on the Sp-Sp mergers, as all other pair combinations were 
found not to produce the strength of PSB galaxy visible to us in the 
data. 

As the simulations, and thus star formation rates, are only 
computed over a total duration of 3 Gyr, we must assume a star 
formation history and galaxy age for the initial build up of stel- 
lar mass in each progenitor galaxy. We assume that by the start of 
the simulation each galaxy has built both its bulge and disk stel- 
lar mass through an exponentially declining star formation rate: 
SFR= SFR +nexp(-t/r) where r = lGyr. SFR is the initial 
SFR of the simulation i.e. ~1.5 Gyr before the merger takes place. 
n is defined such that the galaxies are 5 Gyr old at the start of the 
simulation, i.e. the stellar populations are approximately 7.5 Gyr 
old at the onset of the merger. Both progenitors in a simulation are 
assumed to undergo the same initial SFH. We note that any simi- 
lar star formation history can be assumed with which to build the 
pre-merger galaxies, without significantly affecting their evolution 
in the post-merger post-starburst phase. The creation of the spectra 
from the SFHs of the simulations follows exactly the same method 
as for the toy burst models. 



4.2.2 Post-starburst galaxies in merger simulations 

Figure[6]shows the star formation history and evolutionary tracks in 
PCI/2 for three of the major merger simulations with f gas — 0.2 
(left) and 0.4 (right), BH feedback and prograde, edge-on orbits 
("GO") over the 3 Gyr of the simulation. From top to bottom panel 
the mass of the progenitor galaxies increases. The SFHs of the ma- 
jor merger simulations in general follow similar patterns, with an 
initial weak enhancement in their star formation, quickly followed 
by the starburst and decay to zero star formation. 

We parameterise the simulations using two simple values to 
describe the evolution of the starburst. Firstly, the burst mass frac- 
tion /burst, defined as the mass of stars formed while the SFR is 
at least 1.05 times the continuum rate, divided by the total mass 
of the final merger remnant. The continuum rate is measured in 
the first Gyr of the simulation, and linearly extrapolated across the 
subsequent starburst. The second parameter is the decay time of 
the starburst, tdccay, defined to be the time from the peak of the 
starburst until the SFR is only 10% of the maximum burst height. 
These values are shown in the top left of each SFH panel in Fig- 
ure [6] As found using the toy-models, higher gas fractions lead to 
stronger post-starburst features and longer periods of time spent in 
the post-starburst phase. 

Using this suite of simulations we can investigate several other 
initial conditions which affect the evolution of the stellar popula- 
tion of the merger remnant: progenitor mass ratio; orbital geome- 
try; presence or absence of a black hole. We illustrate the latter two 
in Figure[7] The left two panels show the star formation history and 
evolutionary tracks in PC 1/2 for galaxies with tilted orbital config- 
urations and fg as — 0.2. Comparing to the left panels of Figure 
[6] we see that in this particular suite of simulations, only the pro- 
grade edge-on mergers give rise to post-starburst features. While 
the starbursts are weaker in the other orbital configurations, it is 
likely that the main cause of the difference is in the longer decay 
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Figure 6. Far left: The combined star formation history of both galaxies participating in a major merger during the 3 Gyr SPH simulation. In these panels, the 
galaxies collide in a "GO" orbital configuration and BH feedback is included. Progenitor galaxies have gas fractions of 20% and virial velocities of 80 (top), 
160 (middle) and 320 (lower) km/s corresponding to stellar masses of around 10 10 , 10 11 and 10 12 Mq. The red part of the SFH depicts the "starburst" phase, 
defined to be where the SFR is 5% greater than the continuum SFR. The values in the top left of the panels are the mass fraction of stars formed during the 
starburst and the decay timescale of the starburst (see text). Next left: the corresponding evolutionary tracks in PC1/PC2, with the red section corresponding to 
the red section in the left panel. The dotted box shows the PSB region defined in this paper. The time that the merged galaxy spends in the PSB phase is given 
in the bottom left of the panels. Right pair of panels: The same as the left panels, but for progenitor galaxies with gas fractions of 40%. 




Figure 7. Same as for Figure[6]but for different initial conditions. Left pair of panels: Simulations in which the progenitor galaxies collide with inclinations 
other than edge-on. From top to bottom the orbits are "G7", "G10" and "G13". The galaxies have gas mass fractions of 20% (compare to left hand panels of 
Figure[6] Right pair of panels: Simulations in which the black hole feedback has been switched off. The progenitor galaxies have gas mass fractions of 40% 
(compare to right hand panels of Figure[6). 



times of the star formation rate. However, we caution that the tilted 
orbital geometries tested here have typical tilts of > 30 degrees and 
are retrograde. Other relevant parameters, such as impact parame- 
ter, galactic structure and variations in stellar feedback prescription 
have not been investigated. While the simulations indicate that a 
relatively strong interaction is re quired to produce a post-starburst 
galaxy (see also lCox et"aTT 2008, JNB09), a full study of the post- 
merger stellar populations as a function of impact geometries is 
beyond the scope of this paper. 

The right two panels of Figure [7] show the evolution of the 
stellar population with BH feedback "switched-off" and f gaa = 
0.4. Comparing to the right-hand panels of Figure[6]we see that the 
burst mass fractions and time spent in the post-starburst phase are 



similar with or without the presence of BH feedback. The primary 
causes of the shut-down in star formation are supernovae feedback 
and exhaustion of gas supplies, rather than any additional energy 
supplied by the AGN. In two out of three cases the post-starburst 
features are stronger when BH feedback aids in expelling the gas, 
but the main effect is focused on the later stages of the burst where 
any residual star formation is halted more quickly. Therefore, we 
conclude that from these spectral indices alone, there is no way to 
confirm or refute the suggestion that strong BH feedback provides 
the mechanism for shutting off star formation in the Universe. 

In Figure [8] we summarise the results by showing the time 
spent in the post-starburst phase and maximum value of PC2 
reached after the merger has occurred, for all of the Sp-Sp merger 
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Figure 8. Same as Figure [5] except for SPH galaxy merger simulations with "GO" orbits. Different shapes represent different gas fractions of the progenitor 
galaxies, filled symbols are 1:1 mass ratio mergers and open symbols are 3:1 mass ratio mergers. Increasing symbol size denotes increasing stellar mass. Gray 
symbols are mergers with black hole feedback "switched-off", 40% gas fractions and i> v ; r = 80, 160 and 320km/s. Left: The maximum PC2 attained by the 
models after the starburst has occurred, as a function of burst mass fraction. The dotted line indicates our observational limit above which PSBs are detected. 
Right: The time spent in the PSB phase. 



simulations with prograde in-plane orbits ("GO", i.e. directly edge 
on). The results paint a very similar picture to that shown by 
the toy-models, but we can now associate burst mass fractions to 
physically meaningful parameters such as gas fractions and galaxy 
masses. We find that mergers must induce burst mass fractions of 
at least 5% for them to be detectable in our observationally de- 
fined PSB phase, and the time spent there is greater for more gas 
rich and more massive mergers. For a 3:1 merger (open symbols) 
to appear in the PSB phase, a very high gas fraction of 80% is ap- 
parently required. The duration of the post-starburst phase tends to 
a maximum of 0.6 Gyr, with lower gas fractions leading to shorter 
timescales of 0.3-0.4 Gyr and burst mass fractions of 5-10%. Simi- 
lar to the toy models, we find that decay times of the starburst must 
be shorter than ~0.3 Gyr for a post-starburst galaxy to be observed. 
The results of three simulations with BH feedback "switched-off", 
f 3 as = 0.4 and a range of virial masses, are shown as grey sym- 
bols. These can be compared directly to the filled triangles in the 
same figure, in which BH feedback has occurred. We note that the 
burst mass fractions and starburst decay timescales derived in this 
section a re the same as those fo und for luminous infra-red galaxies 
(LIRGS, iMarcillac et aT]|2006l) . 

Many of the simulated galaxy mergers in the full library of 
79 simulations do not result in a post-starburst remnant according 
to our demanding observational criteria. Even for gas rich major 
mergers, the starbursts can occasionally be weak and/or have long 
decay times, thus not creating the sharp discontinuity in SFR re- 
quired for them to appear in the post-starburst phase. Such galaxies 
may enter the red sequence through the "green valley" if their star 
formation rate continues to decay. 



5 MASS DENSITIES AND MASS FLUXES OF 
POST-STARBURST GALAXffiS AT Z~0.7 

In the preceding sections we have culled a sample of post-starburst 
galaxies with 0.5 < z < 1.0 from the VVDS redshift survey. 
Through comparisons with toy models and SPH galaxy merger 
simulations, we have shown how the spectroscopic features in- 
dicate that the galaxies have undergone a recent major, gas rich 
merger. In this section, we will place the post-starburst galaxies 
into a global context, in terms of their stellar masses, star formation 



rates and their number and mass densities. Finally, we will address 
the interesting question: how much mass enters the red sequence 
through the post-starburst phase? 

5.1 Photometrically derived physical properties 

We can now ask how the mass of the post-starburst galaxies com- 
pares to that of the blue- and red-sequence galaxies? Figure[9]shows 
the mass distributions of our different classes of galaxies, before 
and after corrections for sample selection (TSR, QSR and SSR) 
and survey volume. The galaxies have been classified according to 
their spectral types i.e. position in PC1/PC2. For clarity, we have 
combined the starburst and starforming galaxies into a single class; 
their mass distributions are very similar. The left hand panel shows 
the mass distribution before correcting for completeness effects. 
The right hand panel shows the completeness corrected mass dis- 
tribution, indicating that the post-starburst galaxies primarily have 
log(M/M Q )< 10.5, i.e. to the lower end of galaxies in the red 
sequence and with a similar distribution to galaxies on the blue se- 
quence. It is also worth noting that the completeness corrected mass 
distribution of the post-starburst galaxies continues to rise to lower 
masses, and does not turn over before we reach the completeness 
limit of the sample. Deeper spectroscopic surveys will be required 
to fully probe the mass distribution of post-starburst galaxies. 

Our spectroscopic analysis has focused on uncovering the spe- 
cial class of post-starburst galaxies, and not on recovering further 
physical parameters. As described in Section |2T2l and more fully in 
W08, the multiwavelength (UV-IR) SEDs of these galaxies have 
been analysed to obtain more physical properties than just stel- 
lar mass. In Figure [10] we present the light- weighted mean stel- 
lar age, specific star formation rate (SFR/M*), time of last burst 
and dust content of our full sample of galaxies, classified by their 
PC1/PC2 values. Overall, the relative agreement between the SED 
fitting and VVDS spectral fitting is remarkable: galaxies which are 
spectroscopically classified as starforming are found, from multi- 
wavelength photometry, to be younger, have high star formation 
rates, more recent starbursts and higher dust contents compared 
to galaxies spectroscopically classified as quiescent. The multi- 
wavelength photometry is able to distinguish that the post-starburst 
galaxies have ages and star formation rates intermediate between 
the quiescent and star-forming galaxies. Additionally, the estimated 
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Figure 9. The mass distribution of galaxies, classified into starforming (blue), quiescent (red) and post-starburst (orange, dot-filled) by their spectral type 
i.e. their PC1/PC2 values shown in Figure [3] The orange line-filled histograms indicate the subset of post-starburst galaxies with no ongoing residual star 
formation. The verticle dotted line indicates our post-starburst mass completeness limit. Left: the raw data, with no correction for survey incompleteness 
effects. Right: for galaxies with stellar masses above our completeness limit, after correction for survey incompleteness effects. 



time of their last starburst has a lower mean value than quiescent 
galaxies, peaking at 10 8 -10 9 years, exactly as expected from the 
WDS spectral analysis. We can see that the post-starburst galaxies, 
identified purely through their stellar continuum as having under- 
gone a recent starburst, have a range of specific star formation rates 
which place them at the lower end of the star-forming sequence 
through to the quiescent galaxies. In Section [531 we will use this 
result to set a lower limit on the mass flux onto the red sequence 
through the PSB phase, by identifying those galaxies which have 
completely ceased star formation. 

5.2 Mass and number density 

Summing up all the galaxies in our sample above the mass com- 
pleteness limit of logM = 9.75 with the appropriate weightings, 
we obtain a total stellar mass density of log p* = 8.23^q'qu • The 
errors account for individual errors on the stellar masses and Pois- 
son errors from the SSR and QSR weights. As we have placed 
a stringent criterion on the minimum signal-to-noise ratio of the 
spectra, it is important to ensure that w e recover the correct total 
stellar mass density found in other work. IPozzetti et alj ( 120071) cal- 
culate the total mass density of galaxies as a funct ion of redshift in 
the WDS spectroscopic sample, using the same Chabrier ( 2003) 
IMF as used in the present work. Integrating a Schechter function 
with the parameters[j given by Pozzetti et al. for the redshift range 
0.7 < z < 0.9 above our mass limit of 9.75 gives a stellar mass 
density of log p* — 8.28. Although the statistical errors are always 
< 10%, Pozzetti et al. quote a typical scatter between different 
methods of mass determination and luminosity function estimates 
of at least ~ O.ldex. Taken together with the slightly different red- 
shift range, our results appear to be entirely consistent with those 
of Pozzetti et al.. 

Table [TJ gives the mass and number densities for each of our 
galaxy classes. Note that the completeness of our survey is lower 
for quiescent galaxies than for other classes. 

7 We compare to the /-band selected sample with masses determined from 
complex star formation histories (row 7 of their Table 2), for better compar- 
ison with our own sample and masses 



Focussing on the post-starburst galaxies, we find a num- 
ber density of of 1.0 x 10~ 4 per Mpc 3 and mass density of 
2.69 x 1O 6 M0/Mpc 3 . It is instructive to compare their number 
density to other objects considered to be unusual in the Universe 
at these redshifts. The number density of our post-starburst sam- 
ple is 0.3-0.5dex above the number density of powerful X-ray se- 
lected AGN with 2 -8 KeV luminositie s above 10 43 erg/s at redshifts 
0.5 < z < 1.5 dBargeretalj [2005). They are are factor of 10 
more numerous than star-forming galaxies at these redshifts with 
radio lumin osities similar to local ultra luminous infra-red galaxies 
(ULIRGS. ICowie et alj|2004h . The post-starburst galaxies are con- 
siderably more common than sub-mm selected galaxies with FIR 
luminosities > 6 x 10 Lq at z ~ 0.9, which have n umber densi- 
ties of 3 x 10~ 6 /Mpc 3 /decade JChapman et~alll2005l) . These com- 
parisons are designed for orientation only, and not to imply evolu- 
tionary sequences: the very different duty cycles and "on" times for 
different classes of rare objects make inferences about the relation 
between them difficult. We will return to a quantitative comparison 
with QSO number densities and major merger rates in Section|7] 



5.3 The mass flux of post-starburst galaxies onto the red 
sequence 

We are interested in how much mass may enter the red sequence 
after a starburst and subsequent fast quenching of the star formation 
as seen in simulations of major mergers. We estimate the total mass 
flux which passes through the post-starburst phase to be: 

NpSB , „ 

PA^Q,PSB = > — (3) 

IPSB Wi.Vmax.i 
i=l 

where M* is the stellar mass of the galaxy in solar masses, w is the 
combination of the selection functions w SSK , to TSH and io c * SR ', 
and Vmax is the correction for the volume in which each galaxy 
can be seen. The sum is over post-starburst galaxies above the mass 
completeness limit of log(M/M„)> 9.75. tpsB is the time that a 
galaxy will be seen in the post-starburst phase. From Section|4]we 
find tpsB to be less than 0.6 Gyr. Galaxy mergers with lower gas 
fractions of around 20-40% result in slightly smaller burst mass 
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Figure 10. The physical proper ties of our spectroscopic galaxy sample above a mass limit of log(Af/MQ) > 9.75, derived from fitting their multiwavelength 
SEDs (see Walcher et al. 2008). From top left, clockwise: the light-weighted mean stellar age; the specific star formation rate (SFR/M*); the time of the last 
starburst; the extinction in the V band. The galaxies have been classified by their PCI/2 values into quiescent (red), starforming (cyan), starburst (blue) and 
post-starburst (orange). The histograms represent the raw data, they are not corrected for incompleteness effects. 



Table 1. The mass and number densities of each of our classes of galaxies, for a mass limit of log(M/M ( ,)> 9.75. 



Class 


Quiescent 


Star-forming 


Starburst 


Post-starburst 


Total 


log(Number/Mpc 3 ) 
log(p*/M Q /Mpc 3 ) 


-2.5 

7 so+O-Ol 
'•° y -o.oi 


-2.9 

7 ao+0.02 
'• 8S -0.02 


-3.1 

7 17 +0.03 
' - 1 ' -0.04 


-4.0 

6 43+O.O6 


-2.3 

O 9 n+0.01 

°- zo — 0.01 



fractions and shorter time spent in the post-starburst phase of ~ 
0.35 Gyr. 

During the selection of our post-starburst galaxy sample, we 
have imposed no restrictions on nebular emission line strengths (see 
Section^. However, galaxies which show post-starburst signatures 
as well as nebular emission lines may be in the process of regener- 
ating their star formation and may therefore subsequently return to 
the blue-sequence. In order to derive a firm lower limit on the mass 
flux entering the red sequence through the PSB phase, we elect to 
sum only those post-starburst galaxies which show no evidence of 
star formation. Inspecting the distribution of SSFRs of our post- 
starburst galaxies (top right panel of Figure 1 101) shows the majority 
have SSFRs similar to the lower end of the star-forming sequence, 
6 objects have SSFRs below 10 -11 per yr, 5 of which lie above our 
mass completeness limit. We note that these galaxies also have no 
measurable [Oil] emission. The mass distribution of these galax- 
ies is indicated in Figure [9] by the filled-line histogram and their 
positions marked in PC 1/2 by the open circles in Figure[3] 

Summing Equation [3] over the 5 post-starburst galaxies above 
the mass completeness limit with negligible ongoing star forma- 
tion as derived from their SEDs, assuming an upper limit for tpsB 
of 0.6 Gyr, results in a lower limit on the mass flux entering the 



red sequence through the post-starburst phase of pa^q,psb > 
0.0022l°;°°° 2 M /Mpc 3 /yr. Our best estimate for the mass flux 
through the post-starburst phase onto the red sequence, assuming 
tpsB = 0.35Gyr, is pa^q,psb = 0.00381° ZT M /Mpc 3 /yr. 
We note that this is still a lower limit, as we have excluded all those 
post-starburst galaxies with detectable levels of ongoing star for- 
mation. This residual star formation may still decay as the post- 
starburst ages, we note that the younger post-starburst galaxies 
(those with smaller PCI) are found to be more likely to have resid- 
ual star formation than older post-starburst galaxies (see circled 
points in Figure [3}- If all the post-starburst galaxies were to en- 
ter the red sequence, for tpsB = 0.35Gyr the mass flux would be 
Pa^q,psb = 0.0077 M /Mpc 3 /yr 

To determine how important this mass flux is i n terms of the 
global build up of the red sequence, we compare to lArnouts et al.l 
j2007l) . They use the multiband photometric data in the same 
VVDS-02h field to measure the mass build up of the red sequence; 
comparison with this result thus limits our exposure to cosmic vari- 
ance. The galaxy masses are measured using the same code and 
model dataset as in this paper, thus eliminating one more potential 
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source of systematic error. They finqf|a total mass growth of the red 
sequence of pA—,Q = 0.0098 M0/Mpc 3 /yr. If we assume that all 
galaxies that enter the red sequence subsequently remain on the red 
sequence, then our results show that > 22^g% and likely as much 
as 38tn % of the growth of the red sequence at z < 1 takes place 
from galaxies which have passed through the strong post-starburst 
phase. If all our PSB galaxies with log(M/M Q )> 9.75 enter the 
red sequence, they account for ~ 80% of the growth of the red 
sequence. We note that these numbers are significant, and perhaps 
even surprisingly high. We will return to this point in Section|7] 

The quoted errors on our mass flux include errors on the indi- 
vidual galaxy stellar masses, calculated from the probability dis- 
tribution function resulting from the Bayesian SED fit, and the 
Poisson errors on t he QSR and SSR weights. As described in 
IPozzetti et al.l J2007h . systematic errors caused by the methods used 
to estimate the masses and selection of samples lead to a typi- 
cal scatter in stellar mass densities of 0. 1-0.2 dex. If we discard 
our statistical errors, and assume instead a fixed 0.15 dex error 
on the stellar mass density of the post-starburst galaxies, we find 
Pa^q,psb = 0.0038 tn M w /Mpc 3 /yr However, our compar- 
ison with the results of lArnouts et alj (12007) should be more robust 
due to the use of the same stellar masses, similar survey selection 
criteria and survey area. 



6 COMPARISON AT LOW REDSHIFT 

To create a comparison sample at low redshift, we combine the 
Sloan Digital Sky Survey (SPSS) spec troscopic data release 5 
(DR5. lAdelman-McCarthv & et alj2007» with the U KIRT Infrared 
Deep Sky Survey (UKIDSS. lLawrence et alj|2007l) large area sur- 
vey (LAS), data r elease 3 (Warren et al. in prep.). The SDSS DR5 
galaxy catalogue dStrauss et alj|2002h covers 5740 sq. degrees and 
contains more than 670,000 spectra of galaxies with a median red- 
shift of z ~ 0.1 and tab < 17.77. The accompanying SDSS pho- 
tometric survey provides optical ugriz photometry for each object. 
UKIDSS is an ongoing survey usi ng the UKIRT Wide Field Cam- 
era (WFCAM. ICasali et al.l [2007h to obtain YJHK near-infrared 
photometry covering 4000 sq. degrees in the same region of the 
sky as the SDSS su rvey. The photometric system is described in 
dHewett eTa l. 2006), the pipeline pr ocessing and science archive 
are described in dHamblv et ai]|2008h . The LAS reaches a depth in 
K of 18.4. 

For the purposes of this paper, the SDSS provides the spectro- 
scopic data required to locate the post-starburst galaxies, and the 
near-infrared photometry of the UKIDSS survey is combined with 
the optical photometry of SDSS to obtain accurate stellar masses. 
Requiring the objects to be detected in all YJHK bands gives a 
combined survey area of 916 sq. degrees (S. Maddox, private com- 
munication). 



6.1 Sample selection and incompleteness corrections 

We select SDSS galaxies to have extinction corrected r-band pet- 
rosian magnitudes r\B < 17.7, be primary galaxy targets (PRIM- 
TARGET = TARGETJ3ALAXY or TARGET_GALAXY_BIG or TAR- 
GET_GALAXY_RED), primary catalogue objects (mode = 1) and 



8 They quote a value of pA^.Q = 0.017 ± 0.004 M /Mpc 3 /yr assuming 
a Salpeter IMF, to which they convert their Chabrier IMF stellar masses by 
adding 0.24 dex. Removing this correction results in the value quoted. 



spectroscopically classified as a galaxy (SPECCLASS=2). They are 
matched to UKIDSS galaxies by identifying the closest object 
within 0.5 arcsec. We remove a small fraction of galaxies (3%) with 
bad spectra by imposing a per-pixel-SNR limit of 5, and account for 
this loss in our weighting scheme. Our final sample contains 14822 
galaxies with 0.05 < z < 0.1. 

The redshift range is selected with two things in mind. Firstly, 
we require low enough redshift that our sample is complete down 
to similar masses as for the WDS survey and secondly, we require 
high enough redshift to minimise the aperture bias effect caused by 
the 3" SDSS fibres. This latter property of the SDSS survey means 
that the spectra only probe the central regions of nearby, or mas- 
sive, galaxies, with a strong redshift dependence. For our compar- 
ative census of post-starburst galaxies this is problematic as there 
is currently no known way to robustly detect post-starburst features 
from optical photometry alone. Thus, it is not possible to know for 
sure whether the proceeding starburst was nuclear or global, the 
former being unlikely to be connected to the build up of the red 
sequence. At z — 0.05 the fibre corresponds to a physical diame- 
ter of about 3kpc (h = 0.7), which is small compared to the size 
of a gala xy, but is consider ably larger than the extent of a nuclear 
starburst dBQker et alj |2004). A remaining limitation caused by the 
aperture bias of SDSS, is that post-starburst populations may fail 
to be identified due to population gradients within the galaxies. As- 
suming disks are more likely to undergo starbursts than bulges, this 
would cause us to underestimate the number of PSB galaxies. 

Volume corrections and stellar masses of the galaxies are mea- 
sured in the same way as for the WDS sample (Section l2"3t . using 
the combined UKIDSS and SDSS photometry. One final piece of 
information required is the target sampling rate of SDSS. Although 
the SDSS aims for 100% coverage in the spectroscopic catalogue, 
one of a pair of neighbouring galaxies is occasionally not targeted, 
due to the problem of placing fibres close together on a plate. Sim- 
ilarly, densely populated regions of the sky may suffer from partial 
coverage due to there being insufficient fibres available to cover the 
area. For the purposes of this study, we only require the probabil- 
ity that a galaxy was targeted for spectroscopic follow-up, and no 
further spatial information. We use the SDSS Catalogue Archive 
Server (CAS) to count the number of galaxies within SDSS spec- 
troscopic sectors that are targeted for spectroscopic follow-up, and 
the number with spectroscopic ID numbers (i.e. that were targeted). 
We find that the targeting success rate is 90%. We calculate weights 
for the galaxies by combining this value with the fraction with spec- 
tra above the SNR limit described above. 



6.2 Calculating the Principal Components 

To calculate the Principal Components for the SDSS spectra, we 
first correct for Galactic extinction, convert to air wavelengths, con- 
volve the spectra to WDS resolution and finally rebin onto the 
eigenspectra binning. We then proceed as for the WDS spectra 
(Section l3.lt . to obtain the PCs for the SDSS galaxies. The dis- 
tribution of the first two components is shown in Figure [TT] The 
left panel shows the full sample, and the right panel only galaxies 
with log(M/M Q )> 9.75. Comparing with Figure[3]there are clear 
differences between the SDSS and WDS populations. Firstly the 
tight blue sequence in SDSS contrasts greatly with the cloud in the 
WDS sample. In part this is due to the smaller errors on the spec- 
tral indices, but it may also be due to a more "quiescent" mode of 
star formation in galaxies at low redshifts. The second noticeable 
difference is that the entire population is older (larger PCI). Fi- 
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Figure 11. The first two principal components (PCs) for all galaxies in our SDSS sample (left) and for those galaxies above our mass completeness limit, 
log(M/Mg)> 9.75 (right). For analysis purposes, the samples have been split into quiescent (red), starforming (cyan), star-bursting (blue) and post-starburst 
(orange) classes. These classes are defined in the same way as for the VVDS sample, except for the positioning of the partition between quiescent and active 
galaxies (see Fig.[5J. The median errors of the whole sample (black) and the post-starburst galaxies alone (orange) are indicated in the top left. 
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Figure 12. The spectra of post-starburst galaxies culled from the SDSS 
galaxy sample. The original SDSS spectrum is shown in black, with the 
convolved spectrum in blue (histogram). The red line is the overplotted fit 
of the VVDS derived eigenspectra to the SDSS spectra. 

nally, we can see that there are fewer strong post-starburst galaxies, 
despite the enormous increase in overall sample size. 



6.3 Mass density of post-starburst galaxies at z ~ 0.07 

We identify 6 post-starburst galaxies using the same criteria as 
for the VVDS sample, although from Figure QT] these are clearly 
only the tip of a distinct and identifiable population. The spec- 
tra of the three post-starburst galaxies above the mass limit of 
log(M/M Q )> 9.75 are shown in Figure [T"2l The SDSS sample 
has a high mass completeness even below this mass limit, but for 
the purposes of comparison we retain the same mass limit as for 
the VVDS dataset. The number density of PSBs at z ~ 0.07 is 
5.1 x 10~ 7 , i.e. a factor of 200 lower than in the VVDS survey at 
z ~ 0.7. For comparison, over the same redshift range the number 
density of ultra-l uminous infrared galaxies (ULIRGs) dec reases by 
a factor of ~30 dKim & Sanderslll998l : ICowie et al.ll2004l) . 

The total completeness corrected mass of the SDSS post- 



starburst galaxies with log (M/M @ )> 9.75 is %.&L\% x 1O 1O M , 
giving a mass density of 1.2 ± 0.2 x 10 4 Af Q /Mpc 3 . This is just 
0.43% of the mass density of post-starbursts at z ~ 0.7. 

As in Section |531 an upper limit on the time the post-starbursts 
are visible for is 0.6 Gyr. Following Equation [3] we find a mass 
flux of pa^q,psb > 2.0 ± 0.4 x 10~ 5 M Q /Mpc 3 /yr. Assuming 
our best estimate for the visibility of the post-starburst galaxies of 
0.35 Gyr results in a mass flux of pa->q,psb = 3.4 ± 0.6 x 10 -5 
M /Mpc 3 /yr. 

The SED fitting results indicate that two of the low-redshift 
post-starburst galaxies have no residual star formation (SSFR< 
10~ ri /yr), while one has a SSFR of 10" 10 /yr. However, we cau- 
tion that with only the optical and NIR bands, these values are less 
reliable than for the VVDS galaxies with UV-IR coverage. We have 
therefore used all three galaxies to calculate the present day mass 
flux through the post-starburst phase. The masses of the three SDSS 
post-starburst galaxies are log(M/M Q ) = 10.25, 10.15 and 10.05, 
similar to the masses of the VVDS post-starburst galaxies. 

We can compare our lo w redshift PSB mass flux to that calcu- 
lated by dMartin et al.ll2007t) for galaxies in the green valley. They 
assumed that green valley galaxies defined in NUV-r colours in an 
SDSS/GALEX matched catalogue are all entering the red sequence 
(i.e. all the mass flows from blue-to-red) to calculate a transition 
mass flux of p T = O.O33M /Mpc 3 /yr. Clearly, this value is con- 
siderably greater (a factor of 1,000) than our post-starburst mass 
flux. If the assumption that all green valley galaxies are heading to- 
wards the red-sequence is correct, then the post-starburst pathway 
to the red sequence appears truly unimportant in the present-day 
Universe. Some of this discrepancy may be reduced by appeal- 
ing to the differing m ass limits of the two samples. For example, 
lOuintero et alj d2004l) studied a sample of post-starbur st galaxies in 
the SD SS, using a similar magnitude limit to that of Mar tin et al] 
(2007), which includes many lower mass galaxies than in our sam- 
ple. They measured the "rate density" for post-starburst galaxies in 
the SDSS, i.e. the number of galaxies which pass through the post- 
starburst phase, to be ~ 4 x 10~ 5 /Mpc 3 /Gyr. The equivalent value 
for our high-mass sample is ~ 1.5 x 10 _6 /Mpc 3 /Gyr strongly sug- 
gesting that the measured post-starburst mass flux would increase 
were we to relax our mass limit. 
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7 DISCUSSION: ARE POST-STARBURST GALAXIES 
IMPORTANT? 

In this paper we have presented a sample of galaxies with strong 
Balmer absorption lines and, through comparison with simulations, 
have argued that these are the likely descendants of gas-rich major 
mergers. The strength of their spectral features suggest that they 
are post-starburst galaxies with burst mass fractions of at least 5- 
10%. However, irrespective of whether a starburst has occurred, 
it is certain that their star-formation has been quenched quickly 
with a timescale of <0. 1-0.2 Gyr. As discussed briefly in Section[T] 
such fast quenching times are generally associated with supernovae 
feedback after periods of rapid star formation induced by galaxy 
collisions. An alternative scenario may be rapid fuelin g of star for- 
mation in galaxies from filamentary cold-gas flows ( Ker es et al.l 
120051 : lOcvirketalJ 12008), although direct comparisons with the 
outputs of these simulations would be required to test this. 



7.1 QSOs and post-starburst galaxies 

One popular theory for the origin of the red sequence and the re- 
lation between galaxy bulge mass and central supermassive black 
hole mass is through the triggering of massive outflows driven by 
the central black hole. Such outflows would quickly halt star for- 
mation and therefore the theory suggests an evolutionary link be- 
tween QSOs and post-starburst galaxies. It is therefore of interest 
to compare the number density of the post-starburst galaxies to that 
of QSOs at the same redshift. Any meaningful comparison requires 
knowledge of the QSO lifetime and the visibility time of the post- 
starburst galaxies. The former has b een estimated using a v ariety of 
methods to be tQso ~ 10 7 (see e.g. lGoncalves et al . 2008) and we 
assume tpsB = 0.35 Gyr as previously. Setting 

0qso(M < M x ) _ 0PSB (4) 
£qso tpSB 

we can estimate the magnitude limit of the QSOs, Mx, which 
leads to a space density matching the observed space density of 
post-starburst galaxies of 1 x 10~ 4 /Mpc 3 . Using the parameterised 
QSO luminosity function of ICroomet al. (2004) at z — 0.7, we 
find M x » M*(z = 0.7) + 1.5 = -22.2 where the magnitudes 
are 6^— band and Vega zero-point. That the number density of strong 
post-starburst galaxies coincides with that of moderately powerful 
QSOs, 1.5 magnitudes fainter than M* , does not of course nec- 
essarily imply an evolutionary link. The AGN luminosity function 
continues to rise to fainter magnitudes JHao et al J 1 20050 and un- 
covering the triggering mechanisms of AGN as a function of lu- 
minosity remains one of the biggest observational challenges for 
extragalactic astronomy. 

The number density of the post-starburst galaxies evolves 
rapidl y. Extrapolating the evolution of the QSO luminosity func- 
tion of ICroom eFaT] J2004I) to z = 0.07, and taking iQ S o ~ 10 7 as 
before, we find a QSO number density of 3.1 x 10 -5 Mpc -3 above 
the same magnitude limit of Mx = —22.2, a factor of a few below 
the z — 0.7 number density. So current results suggest that the evo- 
lution in number density of QSOs is not as strong as the evolution 
in the PSB number density. To our knowledge, an observed QSO 
luminosity function at z ~ is not currently available, and would 
be required for a more detailed comparison. 



7.2 Major mergers and post-starburst galaxies 

If, as we have suggested, the strong post-starburst galaxies se- 
lected in this paper are the result of major mergers, it is instructive 
to compare the number of post-starburst galaxies to the number 
of clo se galaxy pairs in the same survey volume. Ide Ravel et al.l 
(2008) measure the fraction of close spectroscopic pairs to de- 
rive the merger rate using the same WDS galaxy sample as in 
this paper. For galaxies with log(M/M„)> 10 they derive a 
merger rate of 2 x 10 _4 Mpc _3 Gyr _1 , with a magnitude differ- 
ence criteria such that they are sensitive to mergers with mass ra- 
tios >4:1. This can be compared with the number of post-starburst 
galaxies, divided by the time during which they are detectable i.e. 
2.9 x 10~ 4 Mpc" 3 Gyr~ 1 for t PSB =0.35 Gyr. Given the small 
number statistics of both samples, uncertainty in both merger and 
PSB timescales and the problem of the unknown magnitude dif- 
ference between the PSB progenitors, it is encouraging that these 
numbers are so close. Of course, larger samples and further detailed 
observations of the post-starburst galaxies will greatly aid our un- 
derstanding of their origin and possible link to gas-rich major merg- 
ers. As with the QSO number density, the current estimates of the 
evolution in merger rate with redshift show l ess evolution than ob- 
served in the PSB number density. Althoug h|de Ravel et ai]( l2008l) 
find that the evolution with redshift is stronger for lower luminosity 
and lower mass galaxies, with the merger rate decreasing by a fac- 
tor of two between a redshift of 1 and 0.5, this is still not as strong 
as the evolution in PSB number density . Finally, we note th at out 
of the 36 close galaxy pairs identified by de Ravel et al. ( 200j|) one 
galaxy is part of the post-starburst sample of this paper. Given the 
small probability of a chance coincidence, this suggests that in at 
least one case the post-starburst stellar population is linked with 
tidal disruption caused by an ongoing major merger. 

The strong evolution in PSB number (and mass) density com- 
pared to the evolution of both QSO number density and galaxy ma- 
jor merger rate, leads us to the question of what may cause such a 
sharp decline in the number of post-starburst galaxies. Through our 
comparison with merger simulations, we have identified two impor- 
tant factors in creating a post-starburst galaxy: gas mass fraction 
and timescale of the starburst. Larger gas mass fractions provide 
more fuel for the starburst, leading to stronger and more prolonged 
post-starburst signatures. Because gas is used-up in the formation 
of stars, it is reasonable to assume that gas mass fraction decreases 
with decreasing redshift, a fact that may play a leading role in the 
global decrease of star formation rate density since z ~ 1. A sec- 
ond effect is the starburst timescale, which is generally assumed to 
be linked to the disk dynamical timescale, which in turn is linked 
to the dynamical time of the halo. This latter value is known to in- 
crease with redshift, by a factor of 2 between a redshift of 1 and 
(E. Neistein priv. comm., iBarkana & LoebfeoOlh . A lengthening 
of the duration of the starbursts, together with a decreasing amount 
of fuel available, may be the dominant mechanism responsible for 
the decrease in number density of post-starburst galaxies. Clearly 
further simulations will be required to test these ideas. 

7.3 The environments of PSB galaxies 

Given the possible link of PSB galaxies to gas-rich major mergers, 
it i s interesting to invest igate the local environments of the galax- 
ies. |Cuctijitietjiy (|200g) measured the density on 5 and 8/i~ 4 Mpc 
scales around galaxies in the WDS survey. Of the 16 PSB galaxies 
with log(M/M Q )> 9.75, 15 have available density measures. A 
control sample of 2122 WDS galaxies with good quality redshift 
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flags was selected for comparison. We do not detect any significant 
difference in the mean or median densities aroun d the PSB galax- 
ies, in agreement with the results o fYanetal.(2008). It is, however, 
interesting to see that the spread in local density values for the PSB 
sample is very large. In detail, 3 lie in underdense environments 
(<5s < 0), 9 in normally overdense environments (0 < 5s < 1, 
while the median value for the control sample is 5s = 0.44 ± 0.02) 
and 3 in strongly overdense environments (<5g > 1). Furthermore, 
there is no clear trend between star-forming vs. non-star-forming 
PSB galaxies. The presence of PSB galaxies in all types of environ- 
ment will be an important constraint for understanding their origin 
and evolution. 

7.4 Building the red sequence through post-starburst 
galaxies 

In Section |531 we compared the mass flux through the post-starburst 
phase to the mass b uild-up of the red sequence as measured by 
lArnouts et al. We found a value of ~40% for PSBs with no 

residual star formation, or 80% for all PSBs with log(M/M )> 
9.75. These numbers, although consistent (i.e. not greater than 
100%), perhaps still appear surprisingly high if we expect other 
quenching mechanisms with slower ti mescales to also play a rol e 
in moving mass onto the red sequence l lvan den Bosch et alj2008l) . 
The key question is whether there is a one-way flow of galaxies 
from the blue to the red sequence? What prevents a galaxy from 
re-starting star formation due to subsequent inflow of gas after ex- 
periencing a major merger? SPH simulations of galaxy mergers in- 
voke strong mechanical AGN feedback to prevent subsequent star 
formation and allow the galaxy to remain on the red-sequence. The 
large amount of mass flowing to the red sequence through the post- 
starburst phase at high redshift, may in fact cause us to question the 
fact that strong mechanical feedback is effective in the long-term 
for the majority of galaxies. Semi-analytic cosmological models 
provide one method to test galaxy formation scenarios in a full 
cosmological context and could help to understand the directions 
of mass fluxes. A full comparison of the recent star formation his- 
tories of the VVDS and SDSS galaxies presented in this paper with 
semi-analytic models is underway. A full census of different types 
of transition galaxies as a function of redshift, when combined with 
accurate measurements of the mass densities on the blue and red se- 
quence, will help to reveal the true importance of mechanical gas 
expulsion mechanisms (feedback) on the global evolution of the 
galaxy population. 

We have found that strong post-starburst galaxies have stellar 
masses similar to those of the least massive galaxies on the red 
sequence. This has important implications for our understanding of 
the physical processes involved in the build-up of the red sequence, 
and in particular the subsequent role of dry-mergers in f orming the 
shape of the red sequen ce mass function observed today dBell et all 
2006j; iNaab et alj|2006l) . 



8 SUMMARY 

Using a PCA analysis of the spectra of the VVDS deep spectro- 
scopic galaxy survey, we have selected a sample of 16 galaxies with 
strong Balmer absorption lines at 0.5 < z < 1.0 and above a mass 
completeness limit of log(M/M Q )> 9.75. Through comparison 
with a suite of SPH merger simulations and toy starburst models 
we have shown that these galaxies are likely to have undergone a 
strong starburst within the last few tenths of a Gyr, with burst mass 



fractions of order 10%, similar to those derived for LIRGS at these 
redshifts (Marc illacetaiT2 006). The key requirement for the obser- 
vation of the post-starburst galaxies, is a fast quenching timescale 
of < 0.1 — 0.2 Gyr. We show that there is a maximum visible life- 
time for the post-starburst galaxies of 0.6 Gyr, but for lower gas 
mass fractions of ~20-40% a lifetime of 0.35 Gyr is more likely. 
The key results of this paper are: 

• Number density: PSB galaxies with log(M/M Q )> 9.75 
have a number density of 10 -4 per Mpc 3 . Assuming they are vis- 
ible for an average of 0.35 Gyr, and that a QSO shines for 10 7 
years, this is equal t o the number dens ity of QSOs brighter than 
M*(z = 0.7) + 1.5 dCroom et al.ll2004f) . 

• Masses and mass density: Summing all the mass in the post- 
starburst population gives amass density of log(p*/M0 /Mpc 3 ) = 
6-431q'q4. The mass distribution of the PSB galaxies rises to the 
mass completeness limit of log(M/M Q )> 9.75, therefore mea- 
suring the complete mass distribution for post-starburst galaxies at 
2 ~ 0.7 will require a deeper survey. The true mass distribution is 
crucial for understanding the relative import ance of different mech- 
anisms for causing post-starburst galaxies dKavirai et alj2007h . 

• Mass flux: We select the 5 PSB galaxies with no residual 
star formation according to multiwavelength SED fitting. Tak- 
ing 0.6 Gyr as an upper limit on the visibility time of the PSB 
galaxies gives a lower limit on the mass flux of pa^q,psb > 
0.0022j;E5 °°° 2 M /Mpc 3 /yr. Taking our best estimate for t PSB of 
0.35 Gyr gives P a->q,psb = 0.0038t°;°°; 4 M a / Mpc 3 /yr. Com- 
parin g this to the rate of build up of the red sequence ( Arno utTet al.l 
2007), we find 38^1! % of the growth of the red sequence at z < 1 
takes place from galaxies which have passed through the strong 
post-starburst phase, assuming all these galaxies subsequently re- 
main on the red sequence. 

• Environment: We u se the local density on 5ft _1 Mpc scales 
from ICucciati et al. to show that PSB galaxies are found in 
all environments, from underdense to strongly overdense. 

• Low vs. High z: We compare our high-redshift results to a 
sample of galaxies with 0.05 < z < 0.1 selected from the SDSS 
and UKIDSS surveys. Despite the enormous increase in sample 
size, only 3 post-starburst galaxies of the same strength as in the 
VVDS sample are detected for the same mass limit. We find the 
mass density of strong post-starburst galaxies decreases by a fac- 
tor of 230 and nu mber density by a factor of 200 from z ~ .7 to 
z ~ 0.07 (see also lHammer et alll997l ; lLe Borgne et alj2006h . 

This paper finds that post-starburst galaxies, although rare in 
the local Universe, are of global importance at higher redshift. With 
larger and deeper spectroscopic surveys, exciting new constraints 
could be obtained on the processes which drive galaxy evolution. 
The joint analysis of imaging and spectroscopy will help to break 
observational degeneracies, especially when combined with "ob- 
servations" of galaxy simulations. On the advent of an era of large 
broad band photometric surveys, it is important not to forget the 
wealth of additional information on galaxy evolution available from 
the additional investment in spectroscopy. 
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